Effects of atomic interactions on Quantum Accelerator Modes. 
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We consider the influence of the inclusion of interatomic interactions on the 5-kicked accelerator 
model. Our analysis concerns in particular quantum accelerator modes, namely quantum ballistic 
transport near quantal resonances. The atomic interaction is modelled by a Gross-Pitaevskii cubic 
nonlinearity, and we address both attractive (focusing) and repulsive (defocusing) cases. The most 
remarkable effect is enhancement or damping of the accelerator modes, depending on the sign of 
the nonlinear parameter. We provide arguments showing that the effect persists beyond mean-field 
description, and lies within the experimentally accessible parameter range. 
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Quantum Accelerator Modes (QAMs) are a manifesta- 
tion of a novel type of quantum ballistic transport (in mo- 
mentum), that has been recently observed in cold atom 
optics In these experiments, ensembles of about 10 7 
cold alkali atoms are cooled in a magnetic-optical trap 
to a temperature of a few microkelvin. After releasing 
the cloud, the atoms are subjected to the joint action of 
the gravity acceleration and a pulsed potential periodic 
in space, generated by a standing electromagnetic wave, 
far-detuned from any atomic transitions. The external 
optical potential is switched on periodically in time and 
the period is much longer than the duration of each pulse. 
For values of the pulse period near to a resonant integer 
multiple of half of a characteristic time Tb (the Talbot 
time [2|), typical of the kind of atoms used, a consider- 
able fraction of the atoms undergo a constant accelera- 
tion with respect to the main cloud, which falls freely 
under gravity and spreads diffusively. 

The non-interacting model is a variant of the well- 
known quantum kicked rotor (KR) [3|, in which the ef- 
fects of a static force, produced by the earth gravitational 
field, are taken into account. The linear potential term 
breaks invariance of the KR hamiltonian under space 
translations. Such an invariance may be recovered by 
moving to a temporal gauge, where momentum is mea- 
sured w.r.t. the free fall: this transformation gets rid of 
the linear term and the new hamiltonian, expressed in 
dimensionless units, reads 
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where p and x are the momentum and position opera- 
tor, k and r are the strength and the temporal period of 
the external kicking potential, g is the gravity accelera- 
tion. The relationship between the rescaled parameters 
and the physical ones, denoted by primes, is k — k'/h, 
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FIG. 1: (Color online) The probability distribution at times 
t = 25 (1 st row) and 45 (2 nd row). ((Red) line: lin- 
ear case (u — 0), (purple) triangles/(green) circles: focus- 
ing/defocusing nonlinearity (u = =pl.25)). In the right col- 
umn enlargements of mode are shown; the position of the 
mode, predicted by ((2)1 is marked by the (blue) vertical dot- 
ted line. 



t = hr'G 2 /M = 4ttt'/T b , rj = Mg'r'/hG and g = rj/r, 
where r\ is the momentum gain over one period, G is 
twice the angular wavenumber of the standing wave of 
the driving potential and M is the mass of the atom. 

Symmetry recovery allows to decompose the 
wavepacket into a bundle of independent rotors 
(whose space coordinate is topologically an angle): this 
Bloch-Wannier fibration plays an important role in the 
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FIG. 3: (Color online) Maximum height reached by the mode 
at time t = 45 as a function of u in position (a) and momen- 
tum (b) representation. In the inset the exponential decrease 
of h m ax for positive u in shown in a semi- logarithmic plot. 



FIG. 2: (Color online) The Husimi function of the QAM at 
time t = 45, in the repulsive (a), linear (b) and attractive case 
(c). 



theory of QAMs [4|. 

QAMs appear when the time gap between kicks ap- 
proaches a principal quantum resonance, i.e. r = 2irl + e, 
with / integer and |e| small. The key theoretical step 
is that in this case the quantum propagator may be 
viewed as the quantization of a classical map, with Id 
playing the role of an effective Planck's constant [J]: 
QAMs are in correspondence with stable periodic orbits 
of such pseudo-classical area-preserving map. We refer 
the reader to the original papers for a full account of the 
theory, we just mention a few remarkable points: stable 
periodic orbits are labelled by their action winding num- 
ber w = j / q, which determines the acceleration of the 
QAM w.r.t. the center of mass distribution 
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The modes are sensitive to the quasimomentum (Bloch 
index induced by spatial periodicity), being enhanced at 
specific, predictable values 0]; also the size of the elliptic 
island around the pseudoclassical stable orbit plays an 
important role (if the size is small compared to |e| the 
mode is not significant [3]). 

We consider in this letter the role of atomic interactions 
in such a system; namely evolution is determined by a 
nonlinear Schrodinger equation with a cubic nonlinearity: 

i<ip(x,t')= \H(t') +u\iP(x,t')\ 2 ] ip{x,t'), (3) 

where u is the rescaled nonlinear parameter, whose sign 
describes an attractive (negative)/repulsive (positive) 
atomic interaction. We will come back to its connec- 
tion with physical units in the end of the paper. The 



condensate wave function is normalized to unity. The 
dynamics does not only acquire in this way a qualitative 
novel form, but, due to the nonlinear term, Bloch de- 
composition into independent rotors breaks down. The 
main scope of this letter will be to numerically scrutinize 
how QAMs are still present in the modified system, and 
explore how nonlinearity modifies their features. In the 
end we will briefly comment upon some stability issues, 
by showing that a more refined description, including loss 
of thermalized particles, does not destroy the scenario we 
get from a mean field description. 

Our analysis will be restricted to QAMs correspond- 
ing to fixed points of period q = 1 of the pseudoclassical 
map; the numerical analysis of nonlinear evolution has 
been performed by using standard time-splitting spec- 
tral methods There are several physical parameters 
characterizing the system: g, r, k and u. Here we mainly 
address the role of nonlinearity u: we fix k = 1.4, I = 1, 
e = — 1, rrj ~ 0.4173, and choose as the initial state 
a symmetric coherent state centered in the stable fixed 
point of the pseudoclassical map (xq ~ 0.3027, po = 0), 
whose corresponding winding number is zero. 

A quite remarkable feature appears when we com- 
pare results for opposite nonlinearity signs (keeping the 
strength |u| fixed), see fig.([I]). As in the linear system, the 
wave packet splits into two well-separated components: 
the accelerator mode (whose acceleration is still compati- 
ble with ([2])) and the remaining part, which moves under 
two competitive contributions, the free fall in the gravi- 
tation field, and the recoil against the accelerating part. 
Note that for the present choice of the parameters, the 
former contribution is negligible compared to the second. 

We remark some features, that are common to what 
we observed for a choice of other parameter values: the 
distribution around the accelerator mode is more peaked 
and narrower in the presence of attractive nonlinearity; 
the opposite happens in the case of a repulsive interac- 
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FIG. 4: (Color online) The probability distribution at t = 45 
for strong nonlinearity (u — ±3); the initial states (see text) 
are shown in the insets (line and symbols as in ng[T}. 
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FIG. 5: (Color online) The quasi-momentum distribution 
function (thick line) for \u\ = 0.5 (a) and \u\ = 3 (b) at 
different times (l)t = 5, (2)t = 25 and (3)t = 45. Pur- 
ple(dark)/green(light) lines refer to attractive/repulsive in- 
teractions. 



tion. This can also be appreciated from a Husimi repre- 
sentation of the modes (see fig.©). 

While for repulsive interactions the spreading of the 
distribution, together with peak damping, seems to de- 
pend monotonically on the nonlinearity strength, the 
attractive case exhibits more complicated features (see 
fig-©)- Enhancement of the accelerator mode is only 
observed for small nonlinearities, while a striking feature 
appears at larger values of \u\, namely the accelerator 
mode is suppressed (see fig.©.)). The intuitive expla- 
nation of this result is that strong focusing nonlinearity 
opposes to the separation of the wave packet into two 
parts; indeed, in the case of exact resonance (namely 
t = 2tt), the mode is absent, so the whole wave freely 
falls without splitting and then the maximum height of 
the wave, plotted vs u as in fig. ((3^) , is found to mono- 
tonically increase to the left towards a saturation value. 
While the behavior shown in fig.© has been observed 
for a variety of other parameter choices, we mention that 
more complex, strongly fluctuating behaviour was some- 
times observed at large focusing nonlinearities. In all 
such cases a bad correspondence between the quantum 
and the pseudoclassical dynamics was also observed, al- 
ready in the linear case. 

We remark that the mode damping is sensitive to the 
choice of the initial state, as shown in fig.©. While 
a gaussian initial wave packet leads to the mentioned 
QAM suppression, we may tailor a QAM enhancing ini- 
tial condition as follows: we take the quasimomcntum 
Po that in the linear case dominates the mode (here 
(3q = tt/t — ti/2 ~ 0.5551 and we drop from the initial 



gaussian all components with \f3 — (3q\ > 0.15. As quasi- 
momentum is the fractional part of momentum, this leads 
to the comb like state of fig.©). Even through quasimo- 
mentum is not conserved due to nonlinearity, the QAM 
is strongly enhanced with respect to the linear case and 
the recoiling part is almost cancelled. 

Another way of looking at the nonlinear evolution with 
techniques that are proper in the linear setting is to con- 
sider the distribution function over quasimomenta, de- 
fined by 



mt) = J2 \{n+p\m)t 



(4) 



This distribution is stationary under linear evolution, its 
shape being determined by the choice of the initial state. 

We consider the evolution of a gaussian wave packet 
(for which the linear / is essentially a constant - the hor- 
izontal red line of fig.©), and probe the effect of nonlin- 
earities of both signs. Typical results are as in fig.©: the 
effect of attractive (repulsive) nonlinearity is to enhance 
(lower) the distribution around a value (3 ~ 0.4. No devi- 
ation occurs for quasimomentum /3q (marked by vertical 
lines), whose wave function, according to fig.©), closely 
follows the linear pseudo-classical island. Again the /3 
peak of the focusing case is suppressed for large focusing 
nonlinearities. 

To make sure that our findings may be experimentally 
significant we discuss some stability issues: the first con- 
cerns decay properties of the QAMs. It is known that 
linear modes decay due to quantum tunnelling out of 
pseudoclassical islands 0]: we checked that, on the avail- 
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FIG. 6: (Color online) (a) Probability inside the island for 
|m|=3 (symbols and line as in fig.flj). (b) The mean number of 
non-condensed particles vs the number of kicks, for u equal to 
0.1, 0.5, 0.75, 1, 2, 5, 7 and 10 (starting from below); 12 terms 
in the sum ([5} are considered. 

able time scale, the nonlinear decay behaves in a similar 
way. In fig.([6H) the probability inside the classical is- 
land is shown as a function of time for the initial state 
of fig-lH})); it has been calculated integrating the Husimi 
distribution of each /3-rotor fiber over the island area and 
summing the contributions of different rotors. 

However in the condensate regime there is another pos- 
sible mechanism that might completely modify the for- 
mer picture, namely depletion of the condensate due to 
proliferation of noncondensed, thermal particles. A stan- 
dard technique to estimate the growth of the number of 
thermal particles is provided by the formalism of Castin 
and Dum which has been employed in similar con- 
texts in |8[. To the lowest order in the perturbation ex- 
pansion and in the limit of zero temperature T — ► 0, the 
number of non-condensed particles is given by: 

(SN(t))=J2M)\v k (t)) (5) 

k 

where Vk(t) is one of the mode functions of the system. 
The modal functions (uk{t),Vk{t)) are pairs of functions 
that represent the time-dependent coefficients of the de- 
composition, in terms of annihilation and creation op- 
erators, of the equation of motion for the field operator 
describing the thermal excitations above the condensate. 
They describe the spatial dependence of these excitations 
and propagate by modified Bogoliubov equations. 

Our findings (see fig.lJSb)) are consistent with a poly- 
nomial growth of noncondensed particles, namely in our 
parameter region (and within the time scale we typically 
consider) no exponential instability takesplace. This is 
consistent with recent experimental work [9( , where 87 Rb 
atom condensate has been used to explore QAMs. In 
a condensate of 50000 Rb atoms with repulsive interac- 
tions is realized. In the case of a "cigar shaped" trap, the 
relationship between the number of atoms in the conden- 



sate N and the effective 1-d nonlinear coupling constant 
u is, in our units, N — ua^_/2ao [lOJ, where ao is the 
3-dimensional scattering length and a± 3> ao is the ra- 
dial extension of the wave function. Using the parameter 
values of the experiment 0], one finds N ~ 10 5 • u and so 
N ~ 50000 corresponds ton~ 0.5. Therefore our range 
of parameters includes the experimental accessible one. 

We have investigated effects of atomic interactions, in 
the form of a cubic nonlincarity, on the problem of quan- 
tum accelerator modes: in particular we have charac- 
terized the consequences of both attractive and repul- 
sive interaction; we have also provided evidences that 
the modes are not strongly unstable when reasonable pa- 
rameters are chosen. 
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